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It is shown in this paper that it is possible, at least in principle, to construct a traversable wormhole 
that is stable to linearized radial perturbations by specifying relatively simple conditions on the 
shape and redshift functions. 
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I. INTRODUCTION 



Wormholes may be defined as handles or tunnels in the 
spacetime topology linking widely separated regions of 
our Universe or of different universes altogether. That 
such wormholes may be traversable by humanoid trav- 
elers was first conjectured by Morris and Thorne [if in 
1988. To hold a wormhole open, violations of certain 
energy conditions must be tolerated. 

Another frequently discussed topic is stability, that 
is, determining whether a wormhole is stable when sub- 
jected to linearized perturbations around a static solu- 
tion. Much of the earlier work concentrated on thin- 
shell Schwarzschild wormholes using the cut-and-paste 
technique [2J. In this paper we are more interested 
in constructing wormhole solutions by matching an in- 
terior traversable wormhole geometry with an exterior 
Schwarzschild vacuum solution and examining the junc- 
tion surface. (For further discussion of this approach, see 
Refs. [3-8].) A linearized stability analysis of thin-shell 
wormholes with a cosmological constant can be found in 
Ref . @ , while Ref . [l0| discusses the stability of phantom 
wormholes. In other, related, studies, the Ellis drainhole 
was found to be unstable to non-linear perturbations [Til ] 
but stable to linear perturbations [12]. According to Refs. 
[3(3) however, such wormholes are actually unstable 
to both types of perturbations. 

A rather different approach to stability analysis is pre- 
sented in Ref. [15j . In Ref. [l6j], an example of a sta- 
ble traversable wormhole connecting two branes in the 
Randall-Sundrum model is considered, while Ref. [17[ 
discusses the instability of scalar wormholes in a cosmo- 
logical setting. 

The purpose of this paper is to show that it is in prin- 
ciple possible to construct a traversable wormhole that 
is stable to linearized radial perturbations. The condi- 
tions on the redshift and shape functions at the junction 
surface are surprisingly simple. 



Using units in which c = G = 1, the interior wormhole 
geometry is given by the following metric [l|: 
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The motivation for this idea is the Schwarzschild line 
element 
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In Eq. |T|), $ = $(r) is referred to as the redshift func- 
tion, which must be everywhere finite to prevent an event 
horizon. The function b — b(r) is usually referred to as 
the shape function. The minimum radius r = r is the 
throat of the wormhole, where b(ro) = Tq. To hold a 
wormhole open, the weak energy condition (WEC) must 
be violated. (The WEC requires the stress-energy ten- 
sor T a p to obey T a ^[i a pfi > for all time-like vectors 
and, by continuity, all null vectors.) As a result, the 
shape function must obey the additional flare-out condi- 
tion b'(ro) < 1 1]. For r > ro, we must have b(r) < r, 
while lim r -j.qo 6(r)/r = (asymptotic flatness). Well 
away from the throat both <!> and b need to be adjusted, 
as we will see. 

The need to violate the WEC was first noted in Ref. [l[ . 
A well-known mechanism for this violation is the Casimir 
effect. Other possibilities are phantom energy [l8[ and 
Chaplygin traversable wormholes (l9j . 

Since the interior wormhole solution is to be matched 
with an exterior Schwarzschild solution at the junction 
surface r = a, denoted by S, our starting point is the 
Darmois-Israel formalism |20l [2l| : if Kij is the extrinsic 
curvature across S (also known as the second fundamen- 
tal form), then the stress-energy tensor S*- is given by 
the Lanczos equations: 



(3) 



where [X] = lim r . 
[Kij] - K± - KT 



. a+ X - lim r ^ a _ X = X+ - X- . So 
ljjj — — ji^, which expresses the discontinuity in 
the second fundamental form, and [K] is the trace of 

[***]■ 

In terms of the energy-density a and the surface pres- 
sure V ', S % a = diag(— cr, V , V). The Lanczos equations now 
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yield 



and 
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A dynamic analysis can be obtained by letting the ra- 
dius r = a be a function of time, as in Ref. 0. According 
to Lobo [ToJ], the components of the extrinsic curvature 
are given by 
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where m s = 47ra 2 <7 is the mass of the junction surface, 
which is a thin shell in Ref. [loj . 

When linearized around a static solution at a = do, 
the solution is stable if, and only if, V(a) has a local 
minimum value of zero at a — ao, that is, V(ao) = and 
V(oo) = 0, and its graph is concave up: V"(ap) > 0. 
For V(a) in Eq. (fl"4"]). these conditions are met [6(. 

Since the junction surface 5* is understood to be well 
away from the throat, we expect a to be positive. 
Eq. (fTU)) then implies that 6(a) < 2M, rather than 
b(a) = 2M, which the Schwarzschild line element ([2]) 
might suggest. (One can also say that the interior and 
exterior regions may be separated by a thin shell. The 
reason for this is that in its most general form, the junc- 
tion formalism joins two distinct spacetime manifolds M + 
and M_ with metrics given in terms of independently de- 
fined coordinate systems x+ and xt Q.) What needs to 
be emphasized is that even if b(a) < 2M, b(a) can be 
arbitrarily close to 2M without affecting the above anal- 
ysis. In particular, V(ao) = and V'(ao) = even if 
lim a _). ao _ 6(a) = 2M, since, by the definition of left-hand 
limit, 6(a) < 2M. The condition V"(a ) > should now 
be written V"(a -) > 0. 
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For future use let us also obtain a': from 
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III. THE LINE ELEMENT 

Given our aim, the construction of a stable wormhole, 
our main requirement can now be stated as follows: apart 
from the usual conditions at the throat, we require that 
6 = 6(r) be an increasing function of r having a contin- 
uous second derivative and reaching a maximum value 
at some r = a. In other words, we require that 6'(r) 
approach zero continuously as r — ¥ a (Fig. 1). Keep- 
ing in mind the Schwarzschild line element we let 
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Again following Lobo [T(| , rewriting Eq. (|10l) in the form 
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will yield the following equation of motion: 




a +V(a) =0. 



(13) 



Here V(a) is the potential, which can be put into the 
following convenient form: 



V{a) = 1 - 



|6(q) + M m 2 ( M — |6(a) 
a 4a 2 



(14) 



FIG. 1: The interior shape function attains a maximum value 
at r = a. 

b(r) = 2M for r > a since M = \b(a). In this manner, 
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both b(r) and b'(r) remain continuous across the junc- 
tion surface r = a. It follows directly from Eq. (fTU)) that 
a = at r = a. It is also desirable to have V = at 
r = a, thereby making S l j — 0. To this end we choose 
$(r) so that 

*'(«-) = 7 ~~7TTT\ - (15) 
v ; o(o-2M) v ; 

[Of course, 5>(r) must still be finite at the throat, while 
$(a-) = $(a+).] For r > a, $(r) = |ln(l - 2M), so 
that $'(a-) = $'(a+). We now have Jf T + - i^V = 
and - if 6 ^ = at r = a. So 7? = at r = a by Eqs. 
as desired. 

For the above choice of $(r), the resulting line element 

is 



So by Eq. ([171). 
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Note especially that 

d . . d . . , d , . d , . 
-r9tt\a-) = —g tt {a+) and — g rr (a-) = — gy r (a+J. 
ar dr dr dr 

(17) 

Since the components of the stress-energy tensor are 
equal to zero at S, the junction is a boundary surface, 
rather than a thin shell [T(| , and K{j is continuous across 
S. 



IV. STABILITY 

As noted at the end of Sec. HH V(ao— ) — and 
V'(oq—) = even if lim a ^ ao _ 6(a) = 2M, since 6(a) < 
2M. In Sec. Mil we saw that in the absence of surface 
stresses our junction is a boundary surface, rather than 
a thin shell: since b'(r) goes to zero continuously as 
r — > a—, b(r) continues smoothly at r = a to become 
2M to the right of a (Fig. 1). This implies that the 
usual thin-shell formalism using the (5-function is not di- 
rectly applicable. To show this, suppose we write the 
derivatives in Eq. (|17p in the following form: 

^9 to, = @(r - a)^g+ v {r) + $[-(r - a)]g~ v (r), 

where O is the Heaviside step function. Then by the 
product rule, 
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+ 6{-(r - a)}— gfiv (a-). 



Up to the second derivatives, then, the (5-function does 
not appear, in agreement with Visser [2lj |: by adopting 
Gaussian normal coordinates, the total stress-energy ten- 
sor, which also depends on the second derivatives of g^ u , 
may be written in the form T^ v — S^S^v + 0(j7)T+ + 
0(— T])T~ , thereby showing the (5-function contribution 
at the location of the thin shell; here a is necessarily 
greater than zero. In our situation, however, Sij — at 
the boundary surface, so that, once again, the (5-function 
does not appear. 

Even more critical in the stability analysis is the need 
to study the second derivative of V(a) in Eq. (fT4")) . Since 
V"(a) involves to" = (iira 2 a)", let us first use Eq. @ to 
write a' in the following form: 



4tt V dr 



.e[-(r-a)]£A* 



As long as b(r) is an increasing function without the as- 
sumed maximum value at r = a, a' will have a jump 
discontinuity at r — a. So a" is equal to 5{r — a) times 
the magnitude of the jump (22j. If 6'( a o) — 0) 011 the 
other hand, the calculations leading to Eq. (fTTjl show 
that a' is continuous at a — clq. It follows that there is 
no 5-function in the expression for V"{clq). 

Without the (5-function, one cannot simply declare 
Aira 2 (j to be the (finite) mass of the spherical surface 
r = a, since the thickness of an ordinary surface is 
undefined. (It is quite another matter to assert that 
dm = Airaa 2 da, which can indeed be integrated over a 
finite interval.) 

Returning to Eq. (JTTJ , when a' is evaluated at the 
static solution, then 6'(&o) = implies that cr'(ao) = 0. 
So a approaches zero, its minimum value, continuously 
as a — > ao~, and, as a consequence, a > in the open in- 
terval (ciq — e, do); here e is arbitrarily small, but finite (as 
opposed to infinitesimal). As a result, a is approximately 
constant, but nonzero, in the boundary layer extending 
from r = eto — e to r — a®. So for a € (ao — e, ao), 
in s = Aifa 2 a is a positive constant, but one that can be 
made as small as we please. Referring back to Eq. (|14l) . 
we now find the second derivative of V, making use of 
the condition b'(a ) = 0. Since m s is fixed, we get 



V"(a -) = - 



b"(a -) b{a -)+2M 



a - 
3m 2 
2(a -Y 



b"{a -)[M-\b{a^)] 



(18) 



Since m s is arbitrarily small, but nonzero, the third term 
on the right-hand side is arbitrarily close to zero, while 
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the last term is equal to zero. From V"(aQ— ) > 0, we 
obtain 



b"(a -) < 



2[b(a -) + 2M] 



(«o-) 2 ' 
Using our arbitrary e, we can also say that 

2[b(a - e) + 2M] 



b"(a ~e) < — 



(«o - e) 2 



The continuity of b(r) and a 2 now implies that 



8M 



b (oq) < o". 



(19) 



This is the stability criterion. 



to an exterior Schwarzschild solution at the junction sur- 
face r — a, where $ and b must meet the conditions 
discussed in Sec. IIIII Our main conclusion is that the 
wormhole is stable to linearized radial perturbations if 
b = b(r) satisfies the following condition at the static so- 
lution a — a$: b"(ao) < —8M/o,q, where M is the total 
mass of the wormhole in one asymptotic region. 

Since the curve b = b(r) is concave down, b"(ao) < 0, 
but its curvature has to be sufficiently large in absolute 
value to overtake 8M/oq — 46(ao)/ag. This condition is 
simple enough to suggest that the form of b(r) can be 
easily adjusted "by hand." 

A function that meets the condition locally can also 
be obtained by converting the above inequality to the 
differential equation 



V. DISCUSSION 

This paper discusses the stability of Morris-Thorne and 
other traversable wormholes, each having the metric 
given by Eq. (TT]), where $(r) and b(r) are the redshift 
and shape functions, respectively. The shape function is 
assumed to satisfy the usual flare-out conditions at the 
throat, while the redshift function is assumed to be fi- 
nite. The interior traversable wormhole solution is joined 



where A is a small positive constant. Confining ourselves 
to the interval (ai,ao], a solution is 



b(r) 



^V^sin ^2^^ k 1 ' 



To the left of a% , b(r) can be joined smoothly to a function 
that meets the required conditions at the throat, thereby 
completing the construction. 
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